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ABSTRACT 



The sky-RMS is the simplest model-independent characterization of a cosmological 
anisotropy signal. The RMS temperature fluctuations determined from the COBE- 
DMR four-year sky maps are frequency independent, consistent with the hypothesis 
that they are cosmological in origin, with a typical amplitude at 7° of ~ 35 ± 2 jjK 
and at 10° of ~ 29 ± 1 /xK. A joint analysis of the 7° and 10° 'cross'-RMS derived 
from the data in both Galactic and Ecliptic coordinates is used to determine the 
RMS quadrupole normalization, Qrms-PSi for a scale-invariant Harrison-Zel'dovich 
power law model. The difference in the inferred normalizations either including or 
excluding the quadrupole is a consequence of the sensitivity of the method to the 
Galaxy contaminated observed sky quadrupole. Whilst there are variations depending 
on the data selection, all results are consistent with an inferred Qrms-PS normalization 
of ~ 18 ± 2 iiK. 



Subject headings: cosmic microwave background — cosmology: observations 
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1. INTRODUCTION 

Previous analysis of results from the COBE-DMR experiment (Smoot et al. 1992; Bennett 
et al. 1994) has unambiguously demonstrated the existence of cosmological anisotropy in the 
cosmic microwave background (CMB). The observed anisotropy is consistent with that predicted 
by models of structure formation with power law initial fluctuations of gaussian distributed 
amplitudes and random phases (Bennett et al. 1996). 

The sky-RMS signal is the simplest statistical measure of cosmological structure that can be 
computed from the sky maps. Additionally, the observed sky-RMS on a given angular scale is, 
in principle, a convenient quantity to use in the normalization of cosmological models. However, 
Wright et al. (1994) and subsequently Banday et al. (1994) demonstrated at length that it is 
essential to account for both instrument specific details, such as the exact beam response function 
(rather than using a gaussian approximation, for example), and data analysis specific details, such 
as the subtraction of the best-fit monopole and dipole from the maps, which perturbs the inferred 
normalization from that derived from standard analytic formulae. 

In this Letter, we compute the sky-RMS values from the complete four-year COBE-DMR 
data set and use these values to infer the normalization of a Harrison-Zel'dovich, n = 1, power law 
model. 

2. GALACTIC CUT 

Emission from the Galactic plane dominates the observed sky signal at DMR frequencies but 
can not be modelled with sufficient accuracy to enable its subtraction. Instead, we excise those 
pixels near the Galactic plane where the GMB anisotropy is necessarily contaminated by such 
emission. For simplicity, previous DMR analyses (although see Hinshaw et al.; 1995), attempted 
to minimize the contribution from Galactic emission by retaining only those pixels for which 
the Galactic latitude | b | > 20°. However, visual inspection of the DMR maps is suggestive of 
additional Galactic emission from regions lying just outside this cut which are known to exhibit 
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strong emission in the infrared part of the spectrum. Cayon &; Smoot (1995) noted two regions 
associated with known Galactic structure which showed non-CMB frequency dependence and were 
therefore hkely to be Galactic in origin. Such features will inevitably contaminate the measured 
sky anisotropy amplitude, and therefore we have attempted to excise them from the four-year 
analysis. 

Kogut et al. (1996a) have demonstrated a correlation between free- free and dust emission as 
traced by the DIRBE long wavelength (100, 140 and 240 /^m) sky maps at high Galactic latitudes. 
Furthermore, Haslam &; Osborne (1987) noted a strong association in the spatial distributions 
of the emission from the inner Galactic plane as recorded by the IRAS 60 ^m survey and radio 
surveys at 2.7 and 5 GHz. This correlation with the thermal emission suggests that the dust 
responsible for the diffuse IR emission in the plane is mixed with ionized gas from extended low 
density HII regions. On the basis of these arguments, we have used the DIRBE 140 /im sky map 
(after a model of zodiacal emission has been removed - Reach et al., 1995) to trace the strong 
Galactic emission around the Galactic plane. To maintain continuity with earlier treatments, we 
also remove all pixels at latitudes | b | < 20°. We find that the additional regions to be excluded 
from our analysis are the 'flares of obscuration' noted by de Vaucouleurs (1956) in Scorpius and 
Ophiuchus in the northern Galactic hemisphere and Taurus and Orion in the south. The new cut 
leaves 3881 surviving pixels in Galactic coordinates, and 3890 in Ecliptic]^ 

3. SKY SIGNAL ESTIMATION 

The COBE-DMR, four- year sky maps are used in the present investigation. Kogut et al. 
(1996b) have demonstrated that the corresponding systematic error contributions are small, and 
we ignore them in this analysis (the 68% upper limit at 7° is ~ 1.4 /iK). A discussion of the effect 

^"A further extension of the cut in which we also exclude all pixels at | b | < 30° has also been 
computed, leaving 3056 and 3039 pixels in Galactic and Ecliptic coordinates respectively. Lists of 
the surviving pixels for each cut are available from the Banday on request. 
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of residual high latitude Galactic emission is deferred until Section 5. However, we note that since 
synchrotron, dust and free-free Galactic emission have spatial morphologies characterized by a 
steeply falling power spectrum (~ - Kogut et al. 1996a), only structure on large angular 
scales is compromised. As a simple initial test of the extent to which Galactic emission affects our 
results, we perform all computations both including and excluding the observed quadrupole. 

The analysis proceeds as described for the two-year data in Banday et al. (1994) where the 
cross-RMS, <^rmS: between two maps a and b is defined in their eqn (1). The data without 
additional smoothing is referred to as 7° smoothing, since the central lobe of the DMR beam is 
approximately described by a 7° FWHM gaussian. For determining the sky-RMS at 10°, the data 
surviving the Galactic cut are smoothed by a 7° FWHM gaussian kernel with uniform weighting. 

Table 1 shows the observed (i^RMS values determined from independent sky maps at 31.5, 53 
and 90 GHz constructed from the first two years and second two years of COBE-T)MR data. The 
results are consistent with the presence of sky signal in all four years of data. We have computed 
a frequency dependence in the ®rms of the form between pairs of frequencies. Since the three 
DMR frequencies span the expected minimum in the Galactic foreground emission (see Fig. 2 in 
Bennett et al. 1992 and Fig. 3 in Kogut et al. 1996a) then an average spectral index /? computed 
over the three frequencies could appear fiat even with considerable contamination and erroneously 
be ascribed to a small Galactic contribution. The lack of any frequency dependence in the (^rms 
(expressed in thermodynamic temperature units) without any correction for Galactic emissions is 
consistent with the dominant signal being cosmic in original. 



4. LIKELIHOOD ANALYSIS 

In this Letter we restrict our attention to the inference of the Qrms-PS normalization for 
the n = 1 power spectrum model. The analytical form of the probability distribution of the 

^^At 53 and 90 GHz, note that this is also consistent with a small and equivalent Galactic 
contribution at each frequency as suggested by the results in Kogut et al. 1996a; 1996c 
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COBE-DMR data analysis specific iSirms statistic is unmanagable. A Monte Carlo approach was 
adopted to generate the 'Sirms distributions for a grid of Qrms-PS values (with 5000 simulations 
used for each value of Qrms-Ps)- Each simulation generates maps of the sky temperature 
distribution by combining a realization of Harrison-Zel'dovich sky anisotropy filtered through 
the COBE-DMR beam (Wright et al. 1994). The expected variance of the spherical harmonic 
coefficients of the CMB temperature in this case is given by < a^^^ > oc l/[£{£ + 1)] (Abbott Sz 
Wise 1984). Noise realizations are based on the appropriate values of the RMS per observation and 
observation patterns of the specific DMR channels. The effects of noise correlations (Lineweaver 
et al. 1994, Gorski et al. 1996a) are not significant for this work and hence are not included in the 
simulations. The Q^ms-PS-dependent statistical means, variances and covariances of the 7° and 
10° (8)RM5 are derived from these Monte Carlo simulations, and used to construct the gaussian 
approximation to the probability distribution of the ®rms- This, together with the measured 
'SiRMS values, defines the likelihood function C{Qrms-Ps)- 

We consider ®rms values computed between the A and B channels at each of the three 
frequencies, between the 53 and 90 GHz (A+B)/2 or sum (S) maps, and between A and B 
maps made by coadding the corresponding sides at the three frequencies. The inverse-noise- 
variance weights for the 31.5, 53 and 90 GHz maps in thermodynamic temperature are then 
(0.090,0.718,0.192) and (0.062,0.588,0.350) for the A and B channels respectively. 

Table 2 summarizes the four-year ®rms values and the inferred maximum likelihood values 
of Qrms-PS- The small observed differences between the Galactic and Ecliptic values are expected 
since there is some noise rebinning depending on the pixelization scheme. The normalizations 
determined at each frequency or from combinations of frequency agree at the 1 a level. These 
values are completely consistent with previous analyses based on either the first year or first two 
years of data, lying between these earlier normalizations, and arc in excellent agreement with 
values derived by more sensitive techniques using the four-year data (see the accompanying papers 
Gorski et al. 1996b; Hinshaw et al. 1996; Wright et al. 1996). 

Since any likelihood technique is only unbiased asymptotically, we have tested the ®rms 
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technique described above with an additional set of simulations with noise levels appropriate to 
the 53S(8)90S combination. There is evidence for a small, statistically significant bias (uncorrected 
in Table 2) to a lower Qf^fjig— PS- Over a range of input normalizations from 15 to 20 fxK, the bias 
is typically ~ 0.2 /xK. The bias is sensitive to the signal-to-noise ratio, which is why Banday et al. 
1994 saw no evidence for this effect in their analysis of the two-year data 

An additional check on the bias correction can be determined from a one-dimensional version 
of the {Qrms-PS, n) matching technique described in Wright et al. (1996). The cosmic part of a 
simulated map can be considered as a function of Qrms-PS- Thus, the i^rms value of a simulation 
and the resultant inferred quadrupole normalization generated by processing through the standard 
maximum likelihood (ML) machinery, {QML,sim)i is itself a function of Qrms-PS- For any given 
sky realization, we can choose a normalization amplitude - Qmatch, sim ~ which results in the same 
inferred normalization as the real data. For 2500 simulations matched to the observed sky signal, 
the mean Qmatch, sim is a good estimate of Qrms-PS- For the 53S(8)90S Galactic combination, we 
find Qrms-PS = 18.41^:1 fiK {i > 3). 

5. DISCUSSION 

A comparison of results including or excluding the quadrupole shows evidence of a systematic 
shift to a lower Qrms-PS normalization when the quadrupole is included. The effect should 
not be considered as surprising, however. Kogut et al. (1996c) have demonstrated that the 
cosmic quadrupole is counter-aligned to the Galactic quadrupole, and thus the quadrupole as 
observed in sky-maps uncorrected for Galactic contamination is not representative of the CMB 
quadrupole. Whilst the shift is within the range of uncertainties expected from noise and cosmic 
variance, it undoubtedly reflects the partial cancellation of Galactic and cosmological quadrupoles. 
Furthermore, since the RMS estimator (at a given resolution) compresses all of the available 
information in the data to one number, it is very sensitive to the £ = 2 mode. A simple calculation 
(not accounting for noise or sky-coverage) for an n = 1 spectrum shows that ~ 23% of the sky 
variance over the spectral range £ G [2,40] is contributed by the quadrupole. 
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The inferred value for Qrms-PS using the (8)rm5 technique is therefore sensitive to the 
choice of quadrupole amphtude in the sky-maps. We have estimated this bias using a set of 
simulations where the power spectrum normalization has been varied over the range 15 - 20 fiK, 
whilst the realization specific quadrupole amplitude was constrained to lie in the range [4, 28] /xK 
(Kogut et al. 1996c). The sample data sets processed through the maximum likelihood machinery 
yielded biases in the range -1.4 to -2.8 /xK, sufficiently large to correct Qrms-PS and reconcile the 
normalizations . 

If we attempt a more specific correction, the foreground emission can be modelled using the 
DIRBE-14Q nm sky map as a template for free-free and dust emission, and a full-sky radio survey 
at 408 MHz (Haslam et al. 1981) to trace the synchrotron contribution (see Kogut et al. 1996b 
for discussion). After subtraction of the appropriate Galactic contribution, the normalization 
amplitudes including or excluding the quadrupole are in better agreement (Fig. 1): correcting 
the DMR maps for Galactic emission increases the CMB quadrupole. The Galactic correction is 
itself subject to considerable uncertainty, and when taken into account all Qrms-PS values are 
consistent with ~ 18 /xK. 



6. CONCLUSIONS 

In summary, we have used the ®rms statistic derived from the COBE-DMR four-year sky 
maps to infer a normalization Qrms-PS ~ 18 /^K (with a typical significance of ~ 9 a) for the 
amplitude of primordial inhomogeneity in the context of a Harrison-Zel'dovich n = 1 model. The 
sensitivity of the iSirms method to the Galaxy contaminated observed sky quadrupole affords a 
reasonable explanation for the difference in inferred Qrms-PS amplitudes when either including 
or excluding the quadrupole. Normalizations inferred from both Galactic and Ecliptic coordinate 
data sets at all frequencies agree at better than the 1 a level. The COBE-DMR four-year sky 
maps are reliable data sets for the investigation of cosmological models. 

We acknowledge the efforts of those contributing to the COBE DMR. COBE is supported by 
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the Office of Space Sciences of NASA Headquarters. AJB thanks Mark Hindmarsh for hospitality 
at the University of Sussex. 
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Table 1: Observed 0rms in thermodynamic temperature between maps of the first two and second 
two years of COBE-DMR data at 31.5, 53 and 90 GHz. The errors are determined from noise 
simulations for a fixed CMB sky reahzation. The frequency dependence of the data, expressed as 
®RMs{i^) oc is indicated between pairs of frequencies. No statistically significant frequency 
dependence is observed. 
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Table 2: (8irm5 combinations used as input to the likelihood (£) analysis and the inferred 
normalization parameter, Qrms-PS-, for an assumed n = 1 spectrum. 
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FIGURE CAPTIONS 



Fig 1. Likelihood curves for Qrms-PS derived using the 53S^90S ^rms combination at 7° and 
10° assuming n = 1. Thick hues represent the analysis based on the Galactic data, thin lines 
correspond to Ecliptic data analysis. The dashed curves are for the quadrupole included {£ > 2) 
case, solid curves exclude the quadrupole {£ > 3). A consistent shift to lower Qrms-PS is evident 
in the quadrupole included fits. The horizontal error bar shows the range of central Qrms—PS 
values determined from the quadrupole included data after plausible Galaxy corrections have been 
applied. All Qrms-PS values inferred from the COBE-DMR data irrespective of data selection 
are consistent with a normalization of ~ 18 /uK (solid circle) . 



